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Abstract: 

We present the extension of two general algorithms for the treatment of infrared singularities arising 
in electroweak corrections to decay processes at next-to-leading order: the dipole subtraction formalism 
and the one-cutoff slicing method. The former is extended to the case of decay kinematics which has not 
been considered in the literature so far. The latter is generalized to production and decay processes with 
more than two charged particles, where new “surface” terms arise. Arbitrary patterns of massive and 
massless external particles are considered, including the treatment of infrared singularities in dimensional 
or mass regularization. As an application of the two techniques we present the calculation of the next- 
to-leading-order QCD and electroweak corrections to the top-quark decay width including all off-shell 
and decay effects of intermediate W bosons. The result, e.g., represents a building block of a future 
calculation of NLO electroweak effects to off-shell top-quark pair (W+W^bb) production. Moreover, this 
calculation can serve as the first step towards an event generator for top-quark decays at next-to-leading 
order accuracy, which can be used to attach top-quark decays to complicated many-particle top-quark 
processes, such as for tt -I- H or tt -I- jets. 
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1 Introduction 


The calculation of radiative corrections to cross-section predictions at next-to-leading order (NLO) 
saw tremendous progress in the past decades leading to a high degree of automation of calculations. 
This advancement was also driven by the development of general techniques for the treatment of infrared 
(IR) singularities. Such divergences arise in the intermediate steps of the calculation involving massless 
states and cancel in the final result by virtue of the Kinoshita-Lee-Nauenberg (KLN) and factorization 
theorems. The various methods that were proposed in the literature can be broadly classified into two 
types: the phase-space-slicing [I}{5] and the subtraction formalisms [6 -20 . Both approaches allow for 


the analytic cancellation of all IR singularities without jeopardising the flexibility of a fully differential 
numerical phase-space integration. 

In this part of the paper we consider the so-called dipole subtraction formalism, which is one of 
the most established subtraction methods. It was continuously extended since its original formulation 
in Ref. 12 , which was restricted to the computation of NLO QCD radiative corrections involving only 


massless partons, to cover all possible cases that are relevant in scattering processes. In Ref. 16 the dipole 
formalism was developed for the calculation of NLO QED corrections where the emission of photons from 
both massless and massive fermions was considered. The extension of the original algorithm for NLO 
QCD calculations to the case with massive partons was worked out in Refs, 
non-collinear-safe observables in QED corrections was presented in Ref. 


17 18 . The treatment of 


19 where the formalism was 


further extended to cover the various fermion-photon splittings not considered in Ref. 
For the case of decay processes, subtraction terms were constructed in Ref. 


16 


21 for the calculation of 


NLO QCD corrections to the radiative decay of the top-quark. These subtraction terms were based on the 
results of Ref. 22 and constructed to fit into the dipole subtraction formalism. Although the discussion 
in Ref. 21 was restricted to the case of the top-quark decay, the subtraction terms given there should 


in principle be applicable to the decay of any strongly interacting particle, however, with the restriction 
that only massless partons appear in the final state. 

So far, no general algorithm exists for the calculation of NLO QED corrections to decay processes using 
the dipole subtraction formalism. In this paper we fill this gap by explicitly constructing such a process- 
independent subtraction term considering both massless and massive final-state fermions and work out 
its extension to the treatment of non-collinear-safe observables. The further extension of the formalism 
to treat massive partons in NLO QCD corrections to decay processes is then fully straightforward. 

The one-cutoff phase-space slicing (OCS) method, considered in the second part of the paper, was 
initially introduced in Ref. for the calculation of NLO QCD corrections to the multi-jet production cross 
section in e“'"e“ annihilation. As such, its application was limited to the case with only massless coloured 
partons in the final state. It was later extended to allow for incoming partons and massive quarks . 
The OCS method is distinguished by the feature that one single Lorentz-invariant cut is introduced to 
isolate the soft and/or collinear regions in the real-emission phase space. Another frequently used slicing 
method, which is not considered in detail in this paper, is the so-called two-cutoff slicing method [5 23 


As a consequence of introducing two independent cuts on energies and angles, however, this approach is 
not manifestly Lorentz invariant. 

In its original formulation, the OCS method is restricted to amplitudes that involve only two sources 
(particles) of IR singularities, a feature that is naturally supported by colour-ordered QCD amplitudes, 
but neither by sub-leading colour structures in QCD nor by IR singularities originating from 17(1) gauge 
bosons, such as the photon in QED. Furthermore, no results employing mass regularization are documented 
in the literature, although it is a common choice in the calculation of electroweak radiative corrections. 
In this paper we extend the OCS formalism in these two respects, i.e. presenting a formulation employing 
mass regularization and dealing with an arbitrary number of soft or collinear singularities. To this end, 
we employ results from the aforementioned dipole subtraction formalism to obtain an approximation 
of the real-emission matrix elements in the singular regions together with the associated phase-space 
parametrizations suitable for the subsequent analytic integration. Furthermore, we discuss the subtleties 
that arise in the case of processes that involve arbitrary many sources of IR singularities. 


At its design luminosity in Run 2, the LHC will be a top-quark factory, copiously producing top-quarks 
in pairs, as single particles, or in association with other particles. As the heaviest elementary particle 
known to date, the top-quark is a unique window to explore the mechanism of electroweak symmetry 
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breaking and to test the Standard Model (SM) on the one hand, and to probe physics beyond the SM 
on the other. Apart from the top-quark mass, key observables in the top sector comprise the top-quark 
width and all accessible production cross sections, for which precise theoretical predictions are required. 
In these predictions the top-quark width Ft plays a double role: as an observable in its own right and as 
an ingredient in cross-section predictions that take into account top-quark decays. In the latter case. Ft 
enters the Breit-Wigner propagator of the top-quark resonance, so that it should be known to the same 

perturbative order as the calculated cross section. _ 

The corrections to the top-quark decay process t —>■ bW"^ are known at NLO 


24 


25 and next-to- 


next-to-leading order (NNLO) [^-29 in QCD and up to NLO for the electroweak (EW) corrections 30- 
34 . However, all prior calculations either consider on-shell W bosons or only partially account for off- 
shell effects by performing an expansion about the intermediate W resonance. As an application of 
the techniques introduced in this work, we present the fully differential calculation of the QCD and 
EW corrections at NLO to the three-body decay of the top-quark, accounting for off-shell effects of the 
intermediate W boson within the complex-mass scheme 35 -37 and retaining a finite bottom-quark mass. 
We consider both the hadronic and semi-leptonic decay channels and further allow for non-collinear-safe 
photon emission (“bare leptons”) for the latter. The presented calculation of the NLO EW corrections 
to the decays t —>■ bW’*' —>■ ht^vn/hqq'^ thus, provides an ingredient to the (not yet existing) calculation 
of NLO EW corrections to the tt production process pp —>■ bbW'''W“ —>■ 6/ including top-quark decays, 
similar to the respective NLO QCD corrections to the top-quark decay in the NLO QCD calculation to 
pp ^ bbW+W- 


6 / 


38 


41 


On the experimental side, the first direct experimental determination of the top-quark decay width 
was performed by the experiments at the Tevatron collider 42 43 . It, however, still suffered from large 


experimental uncertainties. Recently, the first measurement at the LHC was presented by the CMS exper¬ 
iment with the data collected at -yfs = 8 TeV [44| , significantly improving on the previous measurements, 
with a precision at the 10% level. 

This paper is organized as follows: In Sect. we briefly review the general structure of an NLO 
calculation to decay processes to set up the notation and conventions used throughout this work. Section]^ 
first outlines the dipole subtraction formalism in Sect. 3.1 and subsequently discusses its extension to deal 
with decay kinematics. Massive final-state fermions are treated in Sect. |3.2[ and the special case of light 
fermions is covered in Sect. 13.31 where we also discuss the modifications for non-collinear-safe observables. 
The OCS method is presented in Sect. with a general overview of the formalism in Sect. |4.l| explaining 
in particular the origin and treatment of so-called “surface terms” which arise in the case of processes with 
more than two sources of IR singularities. The various cases that arise from the different combinations of 
massless or massive particles in the initial or final states are covered in Sect. 4.2 The application of the 


two techniques to the top-quark decay at NLO is presented in Sect. where the calculational details are 
summarized in Sect. |5.1| In Sect. |5.2| and |5.3| we present the numerical results for the corrections to the 
decay width and the differential distributions, respectively. The conclusions and an outlook are given in 
Sect. [6l 


2 NLO corrections to decay processes—notation and conventions 

The leading-order (LO) partial decay width of a particle with mass TOq decaying into a specific final- 
state configuration comprising n particles is given by 

where A4o is the lowest-order transition amplitude for the I —^ n decay process with ^ representing the 
possible sum (average) over the final (initial) state degrees of freedom and is the symmetry factor 
for identical particles in the final state. Here we have introduced the abbreviation 

d$„ = d(()(pa;Pi,---,Pn) = (27r)'‘d (2.2) 

for the n-particle phase-space measure, where pi and rrii denote the momentum and mass of the i-th 
final-state particle, respectively, and Pa the momentum of the decaying particle, with p)) = (ma,0) in 
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its centre-of-mass frame. Furthermore, we have used the shorthand notation = {pi ,... ,p„} for the 
associated kinematics given by the set of momenta of the particles. 

At NLO, the decay width additionally receives contributions from radiative corrections, 

pNLO ^ p(0) (2.3) 


with 
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d$ 


1 


n+1 


5, 


{n+l} 




(2.4) 


where 5M. and A4i denote the virtual and real-emission matrix elements, respectively. Here we assume 
that the virtual corrections already include the counterterm contributions, so that all ultraviolet singu¬ 
larities are properly cancelled after renormalization. 

The calculation of higher-order corrections involving the emission of massless particles, in general, 
leads to the occurrence of IR singularities in both the virtual and real corrections which cancel in the 
sum of Eq. (2.4) by virtue of the KLN theorem^ Although the cancellation of IR singularities is well 
understood, in practise, the procedure is non-trivial, since the virtual and real corrections are defined 
on different phase spaces. Moreover, the divergences contained in the real corrections are implicitly 
hidden in the phase-space integration, whereas the singularities in the virtual corrections are made explicit 
using regulators when performing the integration over the loop momentum. To accomplish an analytic 
cancellation of all IR singularities therefore requires to make the singularities of the real corrections explicit 
by means of regulators. On the other hand, it is desirable to retain the flexibility of a fully differential 
numerical approach for the phase-space integration. Both the slicing and subtraction techniques provide 
a prescription to perform a fully differential calculation suitable for a numerical evaluation while having 
the full analytic control over the IR singularities and their cancellation. Both approaches exploit the 
factorization properties of the real-emission amplitude in the various singular limits which we briefly 
review in the following for the case of photon radiation off fermions in order to set up our notation. 


For further details we refer to Ref. 16 which we follow closely. Throughout this work we regularize the 


IR singularities using mass regularization where an infinitesimal photon mass is introduced and the 
masses of the light fermions m/ are retained in the divergent logarithms. A prescription for the conversion 
of the results obtained using mass regularization to the corresponding expressions using dimensional 
regularization is given in the Appendix. 

The singular behaviour of emission amplitudes in the soft-photon limit (k —^ 0) can be described by the 
well-known eikonal approximation where the squared real-emission amplitude, summed over the photon 
polarizations can be written as follows. 






fj' 


pfPf' 


{pfk){ppk) 


\M, 


(2.5) 


where e is the elementary charge. Here the summation over /, f extends over all charged external particles 
(of any spin), and the reduced matrix element A4o is evaluated with the set of momenta that is obtained 
by omitting the momentum k of the soft photon. We have further introduced the sign factors af = ±I 
describing the charge flow, which are defined as cr/ = -1-1 for incoming particles and outgoing anti-particles 
and (7/ = — 1 for incoming anti-particles and outgoing particles. Overall charge conservation can then be 
expressed as 

E^/Q/ = 0. (2.6) 

/ 


^In case of processes with identified hadrons in the initial and/or final state, collinear singularities arise for which 
the inclusiveness assumptions of the KLN theorem are not fulfilled. These universal singularities can be absorbed into a 
redefinition of the corresponding (NLO) structure functions given by the parton distribution functions and the fragmentation 
functions, respectively. 
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For the case of decay processes considered here, the mass nia of the decaying particle represents 
the maximal scale in the process and therefore will never be considered in the small-mass limit. As a 
consequence, only final-state collinear singularities arise where the radiated photon becomes collinear to 
a light fermion in the final state. The asymptotic behaviour of the squared transition matrix in this limit 
takes the form 


Pjfc —>-0 ’ 

A-, 


(2.7) 


where Ki denotes the helicity of particle i which might change to any other spin value tk^. If i is a spin-^ 
fermion, the spin can only potentially change sign, i.e. there is just a sum over r = ±. In this case, the 
functions are given by 


9-°+\pi,k) = 


(out) 


Pik 


9t-\p^,k) = 


1 + _ rr4 

I - Zi Pik 




2{p^ky 


with the dimensionless variable 




P° 


Pi + fco 


( 2 . 8 a) 

( 2 . 8 b) 

(2.9) 


representing the momentum fraction carried away by the fermion after the collinear splitting / —>■ / -I- 7 . 

( 2 . 10 ) 


For the case of unpolarized fermions, Eq. (2.7|) reduces to 

|2 


Pik—^0 

A., 




Note that Eq. (2.8) does not describe particles with a spin other than i. However, any subtraction 


formalism based on this asymptotics may still be used for particles with another spin (such as a W boson) 
as long those particles are heavy enough to avoid collinear singularities, because the asymptotic behaviour 


(2.5) in the soft-photon limit is spin independent. 


3 The dipole subtraction formalism 

As was pointed out in the introduction, so far, no general formalism exists yet in the literature for 
the calculation of NLO QED corrections to decay processes within the dipole subtraction formalism. We 
present the construction of such a general subtraction term in the following and begin in Sect. [XT] with a 
brief overview of the formalism and of our approach for its generalization to decay processes. The explicit 
form of the corresponding radiator function and its integrated counterpart is given in Sect. |3.2| for the 


case of a massive final-state fermion. Section 3.3 considers the case of light fermions and the extension to 
non-collinear-safe observables. 

The described formalism was already successfully applied in the evaluation of the NLO EW and 
mixed NNLO QCD/EW corrections to Drell-Yan-like W- and Z-boson production in the resonance region 
presented in Refs. [^|^, where the corrections to the gauge-boson decay processes appear as building 
blocks of the full calculation. 

3.1 Overview of the method 

The general idea of the subtraction formalism can be summarized in the following formula, 

y’d$„+i^iMii'=y'd$„+w iMsubN + j (^j[dk]\Msubfy (3.1) 
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where A4i denotes the transition amplitude of the bremsstrahlung process. The auxiliary term |A^sub|^ 
is specifically tailored to act as a local counterterm to the real-emission corrections. The subtraction 
term |A^sub|^ point-wise mimics the real-emission contributions in all its singular limits, i.e. it satisfies 
the following properties, 


\M 


sub 


Y.\Mi 


for fc —>• 0 or pik 0 , 


(3.2) 


where pi denotes the four-momentum associated with a light fermion in the final state and indicates 
the summation over the photon polarization. The explicit asymptotics in the respective limits are given 


in Eqs. (2.5) and (2.7). The property (3.2) ensures that the first term on the r.h.s. of Eq. (3.1) is free 


of IR singularities. As a result, the integration over the full real-emission phase space d<i)„_|_i can be 
performed without regulators which makes it suitable for the evaluation using numerical methods. The 
singularities originally contained in the real-emission corrections are now fully encapsulated in the second 


term of Eq. (3.1) which is constructed in such a way to allow for a direct analytical integration of the 


singular contributions. To this end, the real-emission phase space is factorized into the singular sub-space 
given by the photonic part [dfc] and the remaining phase space integration d<I)„ associated with the non- 
radiative process. The symbol “(8)” indicates that this phase-space factorization may contain additional 
convolutions or summations. Performing the analytical integration over [dfc] results in an expression where 
the singularities are explicitly regularized in some appropriate regularization scheme. These divergences 
can then be analytically combined with the respective singularities appearing ii^the virtual corrections, 
so that possible cancellations can be made explicit. The remaining integration d$„ is free of singularities 
and can be evaluated numerically. 

The dipole subtraction formalism exploits the factorization of |A^i|^ in the singular limits by 
partitioning the subtraction term |A^sub|^ into so-called dipoles |A^sub,//'|^ according to 

|Afsub(d’r!,-|-l)l^ = ^ |A^sub.//'(d>n+l)l^ 

fj' 

f¥^f' 




fJ' 

f¥^f' 


(3.3) 


where / and f are called the emitter and spectator particles, respectively. The dipole terms are written 
as a product of the auxiliary functions and the squared matrix element |A^o|^! where the former 

represent universal functions that encapsulate the singularity structure and can be integrated analytically 
over the singular photon phase space once and for all. The associated squared matrix element |A^o|^ 
is evaluated on the reduced n-particle phase space ^njf that remains after isolating the singular one- 
particle sub-space [dfc] of the photon. These dipole mappings are required to reproduce the appropriate 
kinematic configuration in the singular regions and further respect exact momentum conservation as well 
as all mass-shell conditions. 

A plain extension in the spirit of the dipole subtraction formalism for the decay process would construct 
two new subtraction terms: a dipole where the decaying particle a is the emitter, and second where 

a acts as the spectator . However, since there are no collinear singularities associated with the initial- 
state particle a, we combine the two subtraction terms to a single term called which 

is illustrated in Fig. H in terms of the two schematic diagrams that would be associated with and 

(sub) 

9a^ ■ 

Such a combined treatment reduces the number of independent dipoles that need to be computed and, 
moreover, allows part of the calculation to be taken over from the scattering case with a final-state emitter 
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P^a = Er=l^ kt 



Figure 1: Schematic diagrams illustrating the subtraction function for the decay of a particle a and 

a charged final-state particle i. The decay dipole merges the two configurations associated with photon 
radiation from the final (left) and initial (right) state into one subtraction term. 


16 


and spectator described in Ref. 
from those of and of 


The asymptotic behaviour of the subtraction function d, 
in the soft limit and explicitly reads 


(sub) 


differs 


16 


ta^l\pz,Pa,k) 

2 piPa mf ml 


fc^O {p^k){pak) {Pik)"^ {Paky 

(3.4) 

i'!^-\pi,Pa,k) 

■ -- 0(1). 

fc-fO 


In the limit of small fermion masses —>■ 0 the asymptotic form in the collinear region is required to be 


.(sub) 

^ia,r 


(Pi,Pa,k) 


Pffc —>-0 


(out) 

9i,r 


iPz,k) 


(3.5) 


with the functions given in Eq. (2.8). 

The subtraction term |Alsub|^ in the dipole subtraction formalism is then given by the following 


expression in the case of decay processes, cf. Eq. ( |3.3[ ), 

|7Wsub(^n+l)|" =5^1-^ subjZa (‘i>n+l)|V5: |-^sub,ij (^n+1) r 




= < aaQadf^^\pi,Pa,k) Mo{^n,ia;TKi) 


+ ajQj k) Mo{‘^n,zj;TKi) 


3 


(3.6) 


where the summation over the indices i,j extends over all charged particles in the final state, a denotes 
the decaying particle, and known subtraction functions of the dipoles with a final-state 


emitter-spectator pair given in Ref. 16 


3.2 Subtraction function for massive final-state fermions 

In the following we present the construction of the subtraction function for the decay of a particle 

a of mass Wq with a massive final-state fermion i of mass rrii. The special case of light fermions can be 
obtained by a systematic expansion in rrii and will be discussed in the next section together with the 
extension to non-collinear-safe observables. ^ 

In the construction of the phase-space mapping ^n,ia for the decay dipole it is convenient to keep the 
four-momentum of the decaying particle a fixed. This constraint can be maintained by the deformation 
of the momenta ki of all particles that are not involved in the respective dipole term via a Lorentz 
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transformation. Such a transformation is necessary in order to compensate the recoil induced by the 
radiated photon so that four-momentum conservation and all mass-shell conditions are retained. To this 
end, we define the auxiliary momentum Pia of the “recoiling system”, see Fig. [T] 

n —1 

Pta=Pi-pt-k^ = Y.^f, (3.7) 

Z =1 


which absorbs the recoil of the radiated photon through the transformation 


_ 


V'A((p*-bfc)2,m2,p2^) 


pl^ _ 


PiaPa 


+ ■ 


ml + Pi - mt 


2ml, 


■Pa, 


(3.8) 


where 


A*a = \{m},ml,Pl), 

and X{x,y,z) denotes the Kallen function defined as 

A(a;, y, z) = -\-y"^ + — 2xy — 2xz — 2yz. 


(3.9) 


(3.10) 


We further introduce the quantity 


2 2 rtl 

=m„ - TO,- - Pv, - : 


(3.11) 


where denotes the infinitesimal photon mass that is introduced as a regulator for the soft singularities 
for later convenience. The phase-space mapping of the decay dipole then explicitly reads 


Pa =Pa, 


nP = T)^ — PP 

TO i^a ^ la, 


k^ = A>^,k^, 


(3.12) 


where fc; denotes the momenta of all remaining particles of the decay process, including all neutral particles. 
The Lorentz transformation can be expressed in terms of the momenta PP and PP as follows, 


_ 11 {Pia + Pia)^ {Pia + Pia) 1 / , Pia Pia,u 
V — g v ~ + p2 


(3.13) 


P2 p. p 

This phase-space parametrization is identical to the one used in the scattering case with a final-state 
emitter and a final-state spectator j of Ref. 16 , where particle i plays the role of the emitter and the 
recoiling system the role of the spectator of mass to| = Pi = Pi- In detail, the following substitutions 
have to be made in order to transfer the quantities of Section 4.1 of Ref. 16 to the kinematics considered 
here, 


^2 p2 

^ ^la, 


Pj ^ Pia, 
Introducing the dipole variables 

2p^k 


Pj ^ Pia, Piii 


Via — 


PI ’ 


■Pa, 


^ia — 




PiPlC 


PiPia + Piak' 


we 


construct the subtraction functions dll'll'^ as follows, 


(3.14) 


(3.15) 


d‘fll\pi,Pa,k) = 


1 


{Pik) Ria{yia) 
ml 


1 - Zia{l - yia) 

4 


1 + 


2 to" 


2 1 


- 1 - 2 - - 
J. ^la 


to : 

Pik 


d[l'{^\pi,Pa,k) = 


Pi Ria{yia) [1- Zia{l- yio)]' 
mf (1 - 2ia)^ ria{yia) 


- d’fl^\pi,Pa,k), 


2{piky 


Ria{yia) 


(3.16) 
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where we have defined the auxiliary functions 






naiVia) = 1 - 


2 J^?a( 2 mf 


Vi. 


^ia 1 Via 

Note the similarity of these subtraction terms to the final-final dipoles of Ref. 
that establish the appropriate behaviour in the soft limit (3.4). In detail, the functions defined in 
Eq. (4.4) of Ref. 


16' 


(3.17) 

(3.18) 

up to modifications 


16 


coincide with up to the first term in the second line of Eq. (3.16), which is 

new, and by an additional factor of (1 + 2rrif/P^g) in the eikonal term (first term in the first line). As 
a consequence, many parts of the calculations can be taken over from Ref. [16| , considerably simplifying 
the analytic integration over the photonic part of the phase space. It is easy to verify that the radiator 
function defined in Eq. (3.16) reproduces the correct asymptotic behaviour in the soft and collinear 

limits given in Eqs. (3.4) and (3.5), respectively. 

The analytic integration of the subtraction functions over the photonic part of the phase space gives 
rise to the integrated counterparts 


= StT^ / [dfc(m^,yia,Zia)] d!lX“’{p^,Pa,k). 


(sub) , 


(3.19) 


The integrated decay dipole can then be written as 


2 

d4*n+l I Af sub,2a(41n+l) I — ^ ^iQi ^aQa 


d$„ D. 


(sub) 


(4>„) 


with a = e^/(47r) and 


D. 


mz 


£“+ (^n) = 2 In I ^ 1 + 3 ln(a 3 ) - 2 ln(l - a^) - 2 In 


("'^g - VWaf - 

P?„ 


Pi + 2mj 


In(ai)In 


\aa-2 


+ 2Li2(ui) +4Li2 i ~\ — 




- 4Li2 (-vV(d^) + - In^(ai) - In(ai) - ^ 


Pi J 2 mf AmfPi 

— ^ ln(a 3 ) + —^ In 


A. 



m^^/Wa{PIa + 2 ™?) 


+ {l-al) 


P.4 

la 

1 

2 


[4 arctan(a3 ) — tt] 


9 9 

mfmt 


_ , 2 m|y^ 

PliPi+mi) a,Pi 


and the abbreviations 

Pi + 2 m? - 

ai = 


P^ - 


P^ 


2 m* 


02 = 


P^ + 

Q.n. ' 


03 = 




-^/Wa 


(3.20) 


P.2 ] 


(3.21) 


(3.22) 


3.3 Light fermions and non-collinear-safe observables 

The subtraction functions given in Eq. (3.16) are already defined in such a way to correctly 

reproduce the asymptotic form in the collinear limit as demanded in Eq. (3.5) for the special of case of 





































small fermion masses rrn —>■ 0. The corresponding expressions for the integrated dipoles are obtained by 
expanding Eq. (3.21) for small rrii and omitting mass-suppressed terms, 


Dt':l\^n)=CiPl.m^) + CiPl,ml)+4U2 


In^ 


p 2 


2 In 


rua 


VWa 


-VWa, 


3 In 


rUc 


. Tnt 
3 


TTla j 3^2 ^ 


^ia- 2 ’ 

where we made use of the auxliliary function 


m 






(3.23) 


(3.24) 


which was introduced in Eq. (3.8) of Ref. 16 to describe the universal singular behaviour of the integrated 
dipole functions 

The corresponding result where the IR singularities are regularized using dimensional regularization 


can be obtained from Eq. (3.23) by performing a transition between the regularization schemes using the 


results of Ref. 47 which is briefly summarized in the Appendix. 


In non-collinear-safe observables, photons arbitrarily close to an outgoing charged particle are not 
treated inclusively, but it is assumed that such a configuration can be experimentally resolved. This is, 
for instance, the case for muons in the final state, which are observed in the muon chambers further out 
in the detector, whereas the photons are detected independently in the electromagnetic calorimeter. Such 
an observable definition is sensitive to the details of the collinear / —>■ / -I- 7 splitting. The subtraction 
formalism can be extended to cover this case by exploiting the information on the dipole variable Zia 


which corresponds to the momentum fraction Zi of Eq. (2.9) carried away by the emitter fermion in the 


collinear limit. This is accomplished by treating the n-particle kinematics of the dipole phase space 
as an (n -I- l)-particle event with momenta pi —)■ ZiaPi and fc —)■ (1 — Zia)pi- Eor further details we refer 
to Ref. where this extension was introduced. Owing to the explicit dependence of the observables on 
the dipole variable, the integration over Zia = z. cannot be performed analytically, but is left open for a 
numerical evaluation. Employing a plus prescription, the singular endpoint part can be isolated so that 
the extension to the non-collinear-safe case amounts to an additional term which vanishes for inclusive 
observables, 


d4*n+l I Alsub.ia (d^n-t-1) I ©cut ipi — Z^Q^pi^ k — (1 Zia)pi.^ 


a 


= -- a.Q. aaQa / Cl$„ / dz ^ 7?r;^($„)<5(l - z) + I?):7(4>„, z) 


X 


M 


($„; TKi) ©cut (pi = zp^, k = {l- z)pi, {ki}^ , 


where the plus prescription is defined by 

pi 


fO 


dx[f{x)]+g{x) = ( dx f{x)[g{x) - g{l)\ 
do 


(3.25) 


(3.26) 
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for any smooth test function g{x). The function is obtained after carrying out the analytical 

integration over the dipole variable yia and explicitly reads 


V, 


(sub) 

ia,+ 


{z) = - 


1 + z^ 
1 - z 


In 




+ (1 + z) In 1 + 


\ za 
2 


^[1 - viz)] + (1 + 2) ln(l -z)- 


2z 


K. 


- 2 In 1 - 


viz) 

2zri{z) 

1 + a{z) 


(1 - z)cr(z) 


In 1 + 


1 - z 
7?(z)[l - z?7(z)] 

mL(i - ^) 


t{z) In 




■qiz 


-M-z) 


I 2/rf^(l-2z) ^^/^2/xf„(l-z) + ?7(z)[l + cr(z)] 


Pl<^iz)^ 


r(z) 


‘^Viai^ - z) + qiz)[l - aiz)] 


viz) 


1 - qiz) 


l-z[l-j/+(z)] I inf^ + q{z)) il - z) 




(*) = 1 - 


where y-\-{z) denotes the upper boundary of the jjia integration and is given by 


•n -1 


y+iz)= Ciz) + 1 + + 2 ] 

We have further introduced the dimensionless variable 




2 z(l — z) 


p2 

2 ^ _M 

H'la p2 


and the auxiliary functions 

cr(z) = yjl + ipf^zil - z), 


viz) = 


[^-y+iz)]z for z < i, 

[1 - y+W] (1 - 2=) for i < z. 


(3.27) 


(3.28) 


(3.29) 


(3.30) 


4 The one-cutoff phase-space slicing method 


The main idea of the phase-space slicing approach is the isolation of resolved hard emissions from the 
unresolved soft and/or collinear regions in the real-emission phase space upon introducing a technical cut 
parameter. The hard-emission contributions are free of singularities, so that the corresponding phase- 
space integration can be performed numerically using Monte Carlo methods. The contributions from the 
unresolved regions encapsulate the IR singularities of the real corrections, and the integrations over the 
sub-spaces containing the singularities are performed analytically, exploiting the universal factorization 
properties of amplitudes in the soft or collinear limits. This makes the IR singularities that were implicitly 
contained in the real-emission phase space explicit by means of regulators and allows for an analytic 
cancellation of all IR singularities. Various approaches have been proposed for the introduction of the 
technical cut. The OCS method considered in this work imposes a single cut 2piPj > As on the invariant 
scalar products of the momenta Pi^Pj of two particles that can cause a soft singularity [pi —>■ 0 or pj —>■ 0) or 
a collinear singularity ipiPj —> 0). The method, thus, isolates soft and collinear singularities simultaneously 
and results in a formalism that is manifestly Lorentz invariant. The technical cut parameter As should be 
chosen small enough to suppress all effects of order 0(As) to a numerically irrelevant level. Despite the 
drawback of introducing an arbitrary resolution parameter, such as As, slicing methods offer the possibility 
to widely suppress negative weights and to avoid unbounded weights in Monte Carlo integrations, which 
considerably simplifies a later extension of the Monte Carlo integrator to a Monte Carlo event generator 
with unweighted events. 

Although the OCS method has been extensively studied in NLO QCD calculations involving mass¬ 
less or massive [4 23 48 cases, there is so far no exhaustive survey of results covering all relevant 
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configurations of massive particles in production and decay processes. The OCS method is particularly 
suited to NLO QCD calculations in leading-colour approximation where the corrections are decomposed 
into colour-ordered contributions that naturally involve only singularities associated with the two neigh¬ 
bouring hard radiators of the unresolved parton in the colour ordering. If more than two collinear singu¬ 
larities appear in one contribution, as for instance in the sub-leading colour contributions, the application 
of OCS becomes subtle. In this case, process-specific treatments are required, such as reconstructing the 
IR singularities from the leading-colour building blocks as described in Refs. [^[^. 

In the following, we describe a generic method to apply OCS in the presence of arbitrarily many 
collinear singularities in phase space. Moreover, we are not aware of any application of the OCS method 
to QED radiation processes, where e.g. no “colour ordering” exists leading to amplitude contributions 
with only two collinear singularities. Furthermore, the results in the literature all employ dimensional 
regularization with either vanishing masses or large mass parameters; in QED applications, however, the 
case of small, but finite masses frequently occurs as well (e.g. for electrons or muons emitting photons). It 
should be noted that the case of a small mass m cannot be obtained upon simply taking the limit m —>■ 0 
in the OCS method, since the hierarchy As ^ rri^ is used in the integration over the unresolved regions 
for a radiating particle of mass m. In order to isolate the mass singularity regulated by a small mass 
m (which defines a mass regularization), however, the inverted hierarchy <C As is required. In the 
following we describe OCS for QED radiation covering all relevant mass configurations (with large and 
small masses). The transfer to QCD calculations is straightforward and merely requires the reintroduction 
of colour charges known from the massless QCD casej^ 

4.1 Overview of the method 

The photonic real-emission corrections are divided into the finite “hard” contribution In and the 
singular “soft-collinear” contribution Jsc as follows, 

/d$„+i ^|Mi($„+i)|' = Jh + /sc. ( 4 . 1 ) 

A, 

The two regions are separated with the resolution parameter As imposing conditions Sf^ ^ As on the 
Lorentz-invariant quantities sjj defined as 




(4.2) 


where pf and k are the momenta of a charged fermion / and the photon, respectively. The definition (4.21 
holds irrespective whether the fermion is massless or massive. In order to formally partition the phase 
space into hard and soft-collinear regions, we introduce shorthands for the 0-functions 


9/ = ^(s/7 - ^s). 


% = !“%= “ s/7)- 


(4.3) 


For the hard region we require Sf^ > As for all charged fermions / of the process, so that the 
corresponding contribution to the real correction reads 


la = f d$„+i 


0 = n^/- (4.4) 

/ 


The phase-space integration in In can be performed numerically, producing a result that depends on As 
via terms proportional to In As or In^ As for small As. 

The integral over the soft-collinear region is given by 


/sc = / d<i>„+i ^ |Mi($„+i)|2 (1 - 0). 

A., 


(4.5) 


^In this paper, we consider only photon bremsstrahlung in detail, which up to colour factors literally translates into 
gluon radiation off massless or massive quarks in QCD (with splittings q q*g, g** qg and likewise for q). Other QCD 
splittings, such as g ^ g*g, g ^ 9*9, etc., either are covered by the known results for massless quarks or do not need any 
special treatment for massive quarks, because no soft or collinear singularity exists in those cases. 
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In this integral, at least one of the invariants Sf^ is small, i.e. < As. This means that in the 
limit As —>■ 0, non-vanishing contributions from the soft-collinear region to the cross section originate 
entirely from those terms in the matrix element that lead to soft or collinear singularities, because regular 
contributions tend to zero with As —>■ 0. Thus, in the calculation of soft-collinear contributions to the 
cross section we can employ the universal factorization form of the squared amplitude, as for instance 
described in Section^ In fact it is convenient to directly use |A^subP = l-^subj/'P as defined in 

dipole subtraction in rerms of emitter-spectator pairs. We, thus, can evaluate Isc as follows. 


he 



|Afsub,//'P(l - ©)■ 


For each emitter-spectator pair //' we can write 0 as 


(4.6) 




1 if at least one Sfi^ < As for /" ^ f, /', 
0 otherwise. 


(4.7) 


Using 9f + 9f = 1, we can rewrite the factor (1 — 0) to 

1-Q = {9f + 9f){9f< + 9f,) - 9f9f>{l - Off,) = 9f9f, + 9f9f> + 9f>9f + 9f9f,Qff> 
and decompose he into soft (S), collinear (C), and quasi-soft (QS) contributions as follows, 

he = X! hejf = X! ihjf + lejf + iQsjf), 

///' 


h.ff = J d^n+l\Msnh,ff'\^0fdf', 
lejf = J d^n+i\Msnhjf'\'^{df9f> +9f,9f), 


(4.8) 

(4.9) 

(4.10) 

(4.11) 

(4.12) 


The definition of the various sub-contributions is schematically illustrated in Fig. The integrals hjf 
and lejf are confined to the soft and collinear regions, respectively, and involve only the technical cuts on 
the invariants Sf-y and Sfiy of the emitter and spectator fermions. For this reason, the integrations in hjf 
and lejf over the singular subspaces can be performed analytically for small As, as outlined below. On 
the other hand, the quasi-soft integrals iQsjf involve cuts on all charged particles, rendering a process- 
independent analytical treatment impossible; nevertheless their numerical evaluation is straightforward, 
as outlined at the end of this section. Note that the quasi-soft integrals Iqsjf do not occur if only two 
sources of collinear singularities exist in the bremsstrahlung phase space. It is the occurrence of these 
subtle terms that requires process-specific and non-generic manipulations of the formalism as described 
in Refs. for the sub-leading colour contributions. 

The analytical calculation of the singular phase-space integrals in hjf and lejf starts by using the 
factorized matrix element of Eq. (3.3) and the corresponding factorized (n -I- l)-particle phase space, the 
construction of which is described in Ref. 16 in detail. It is convenient to directly calculate the sum of 
the two contributions. 


Isingjf = hjf + hjf 


= -e^^fQf^fQf J d^njf® J d[k]gf^}'\pf,pf>,k) Mo{<^njf ) 

X {9f9f, + 9f9f, + 9p9f). (4.13) 

In the following we restrict ourselves to the case of unpolarized fermions with the radiator functions 

(4.14) 
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^sing,//' 




Figure 2: Schematic representation of the different contributions to the soft-collinear region for a single 
pair of particles, //'. 


Note that the spin information of a massive decaying particle is still fully accessible in spite of this 
choice, because the introduction of a separate spin-flip function was motivated by the massless 

limit for an emitter fermion, which is discontinuous. For collinear singularities of initial-state fermions, 
the multiplication symbol “(8)” in Eq. (4.13) indicates that a convolution over the collision energy after 
photon emission is involved in the phase-space factorization. In this case, the analytical integration 


over d[fc] can only be performed partially, as described in Sects. 4.2.2 and 4.2.3 below. For final-state 
singularities, however, the squared matrix element |A4o(^’n,//')| is constant in the integration over d[fc], 
so that the singular integration is entirely contained in the universal factor 


Gff'isff';As) = 


d[fc] {9f9f> 


9f9fi + 9f'9f). 


(4.15) 


For initial-state singularities, this factor comprises an essential ingredient as well. Results on Gff for all 
relevant mass configurations are provided in the following sections. 

Finally, we describe a simple way to numerically evaluate the quasi-soft integrals Iqsjp- Recall that 
the integral receives only support outside the soft-collinear region of the emitter-spectator pair //' on a 
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phase-space volume of 0(As), because the photon has to be collinear to at least one other charged fermion. 
Non-vanishing contributions to Iqsjf' in the limit As —>■ 0 can, thus, only result from an enhancement 
in the integrand lifting the As suppression in phase space. Power-counting in the photon momentum in 
the integrand reveals that this can only happen if the photon momentum becomes soft, i.e. the photon 
momentum has to be of order C>(As) in order to lie outside the soft-collinear region of //'. The integral 
Iqsjf is therefore effectively a surface integral of the soft-collinear region of //' and can be evaluated 
upon taking into account only those parts in the integrand that produce soft singularities, i.e. the “eikonal 
terms” in are sufficient. 


(sub.eik) / , \ 1 

9ff' = ^ 


2pfpf, 


71 : 1 % 


Pfk + pfik pfk 

After symmetrizing with respect to the interchange / O /', we can identify the usual eikonal factor 

(sub,eik) 


(4.16) 


(eik) / ; \ (sub.eik) / 7 \ , isuD.eiKi/ 7 \ 

9)f' {Pf,Pf',k) = - g)f, '[Pf,Pf',k) + g):,j: ’ '{pp,pf,k) 


1 r 

2 

PfPf _ 

(pfk){pf'k) 2(p/fc)2 2{pf>ky 


m 


f 


nrii 


(4.17) 


Since only the symmetric combination of g^ff!^\pf,Pf', k) is relevant in the overall quasi-soft contribution 
(summed over all emitter-spectator combinations), Iqsjf can be evaluated via 

Iqsjp =^fQf (^fQf J <^^n+i9ff\pf,Pf,k) Moi^njf) (4.18) 

Note that this form is valid for any combination of incoming and outgoing emitters and spectators. In 
the form given in Eq. (4.181, each term Iqsjf involves the phase-space embedding —>• $«,//' that 

is specific to each emitter-spectator pair. Since, however, the integral receives contributions only from 
soft photon momenta k, we can simplify the evaluation even further. For the relevant contribution we can 
factorize the (n -|- l)-particle phase space into the corresponding n-particle phase space for all particles 
but the photon and the 1-particle photon phase space, leading to the alternative form 


iQSjf =-e'^(TfQfafQf, J d$„ \Mo{<I>n)\‘^ J (2^)32fc0 9ff\pf,Pf,k) 


(4.19) 

which does not make use of any special phase-space mapping for ff' pairs and is, thus, much simpler to 
evaluate than Eq. (4.18). In this approximation, some care has to be taken with respect to the upper 
bound of the photon energy which has to be large enough to cover the full relevant soft region, but 
still small enough to avoid artefacts from terms of the order ©(As In 

In the following section we carry out the integrations over the soft-collinear regions in Aing,//' for the 
various emitter-spectator pairs //' in the initial/final states. 

4.2 Integration over the singular regions 

4-2.1 Final-state emitter and final-state spectator 


Following the notation and conventions of Refs. 16 19 , in the case of a final-state emitter i and a 
final-state spectator j (illustrated in Fig. the full integral over the soft-collinear regions is given by 


G,,(p 2 ,As) = 


p4 

2 yjXij , 


dyzi (1 - y^j) 


r^2{Vi: 




9lT'^\pi,Pj,k) {8^9j -\- 6.^9j + 9,9j), (4.20) 


where 


Pij — Pi Pj 


i ~)‘2 td 2 2 2 2 

Pij Pij ^i ^j ; 






(4.21) 
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Figure 3: Generic diagram with a final-state emitter i and a final-state spectator j. 


with A(a,b,c) defined in Eq. (3.101. The variables 


Pik 


PiPj 


PtPj + Pik + pjk' PiPj+Pjk 

parametrize the photon phase space. Their boundaries are given by 


(4.22) 


Vi = 


2mim^ 

y 


2m i 


?/2 = 1 - 


'Plj - rrij 


P .2 

IJ 


ziAy) = 


(1 - j/)(2m| -H Pf^y) T V- vl^/KjRijjy) 

2 (1 -?/)(to-+to2 -P^2 ) 


where 


R^Jiy) = 


(2m2+i^-/^2/)2-4i^ 






(4.23) 


(4.24) 


Note that the infinitesimal photon mass m^ acts as regulator for the soft singularity appearing at pij —>■ 0. 
The integrand g^^^^\pi^pj, k) is given by 


9ti^'°\puP].k) = 


2 1 


{Pik^Rij {jjij ) 

and the soft and collinear regions are characterized as follows 


2 ^ m:^ 

1 ^ij (1 yij^ Pik 


soft region: 

i-collinear region: 

j-collinear region: 


0 


k ^ 0, k \\ Pi 


k ^ 0, 


Pi 


^27 

®t7 


O 27 

®t7 


^27 

®t7 


< As 


Vij < Ay, 

< As 


Zij > 1 — Az, 

< As 


Vij < Ay, 

> As 


Zij < 1 - Az, 

> As 


Vij > Ay, 

< As 


Zij > 1 — Az, 


where 


Ay = 


As 

W' 


Az = 


Ay 

1 - Vij' 


(4.25) 

(4.26a) 

(4.26b) 

(4.26c) 

(4.27) 


Depending on the masses of the final-state particles (massive or light emitter/spectator), we obtain dif¬ 
ferent results for Gij{P^p As), which are listed below. 


The full integral over the soft-collinear regions, as defined in Eq. (4.20), is sufficient to cover all 
kinematic configurations in which the photon is recombined with collinear fermions in the event selection, 
i.e. in the case of collinear-safe observables. If the photon and the emitter i become collinear, this 

















photon recombination ensures that all *-collinear configurations integrated over in Gij , As) are treated 
inclusively, i.e. either the full contribution represented by Gij{P^j,As) contributes or it is completely 
excluded by cuts. This even holds true for contributions in Gij{Pfp As) where the photon is collinear 
to the spectator j, since the photon has to be soft in this case, i.e. Sjj = 0{As). In non-collinear-safe 
observables, however, the photon can be separated from collinear fermions, which is for instance the case 
for muons detected in the muon chambers. In this case the differential information on the energy flow 
inside the i-collinear cone should be kept in the phase-space integral (4.20). This energy flow in the 
soft-collinear region is controlled by the variable Zij, where it is effectively given by zt = -|- /c°). 

In this situation, the integration over Zij in Eq. (4.20) should be included in the numerical phase-space 


integration, so that instead of the full integral Gy As) the following function oi z = Zij appears 


p.4 


rV2{z) 
2\/^ Jyi(z) 


z) = ^ ^ I dj/y (I - 2/y) fc) {6,9 J 


9i9j OjOi), 


(4.28) 


where the functions yi^ 2 {z) result from inverting the parametrization of the integral boundary defined 
in Eq. (4.23). The evaluation of ^y (P? , As, z) essentially follows the analogous treatment described for 
dipole subtraction in the case of non-collinear observables in Ref. 19 . The procedure is greatly simplified 
upon introducing a plus distribution for the z-integration. 


z) = [eiy (p2, As, z)] -k (5(1 - z)Gy (E’^-, As), 


(4.29) 


which isolates the soft singularity in the endpoint function Gij{Pfj, As), so that [C/y (P?, As, z)] can 
be evaluated for = 0. The explicit evaluation of Qij{Pfj, As, z) and of the full integral Gij{P{^j, As) 
requires a discrimination of the different mass patterns in the emitter-spectator pairs. 


(i) Massive emitter and massive spectator 

For massive emitter and spectator fermions, we have the following hierarchy of parameters, 

<C As <C TO^, m'^j,Pfj. (4.30) 

The full integral Gij{Pfj, As) is tedious, but straightforward. The result involves a discrimination between 
the variables , Pf^ and is given by 


Gy(P)2 ,As)= 


In oi In 


/ m^Ay- \ , , 1 ,„ 2 . 


{As^PiJ 


+ 2Li2(ai)--—h — In 


— 2 In 


As 

m^rrii 


+ hjAPfj) T S{mj+m] - Pfj) f,j^±{Pf^) for rrii^mj. 


with the auxiliary functions 


(4.31) 





(4.32) 


and the variable 


Oi = 


Pjj + - ^/^j 

P2- -k 2mj + 


(4.33) 


To calculate the function Qij{P^j, As, z) for non-collinear-safe observables, we have to invert the func¬ 
tions zi^ 2 {y) of Eq. (4.23) using the hierarchy (4.30). Making use of the fact that only yy values of 0{Ay) 
with Ay <gC 1 are relevant, the functions yi^ 2 {z) are given by 


yiiz) = -r 


2 m 7 


_ 

ij 


i^-z), 


y2{z) = ^ ,^^-^ {l-z). 
Pij ■ 


(4.34) 
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This shows that for mi,mj obeying Eq. (4.30) all z values in the soft-collinear region are of the order 


1 — 0{Ay). The integral As), thus, only receives contributions from a tiny neighbourhood of the 

point {yij,Zij) = ( 0 , 1 ), i.e. it does not involve an integration over the momentum flow in the collinear 
regions. In this case non-collinear-safe observable can be calculated using , As) without the need 

to introduce Qij{Pfj,As, z). 


(ii) Massive emitter i and light spectator j 

Due to the small mass of the spectator, the hierarchy between the parameters is 


^ As ^ m^,P^j, rrij = 0. (4.35) 

The spectator mass can be set to zero, since no collinear singularity exists for the spectator. In fact 
the integral Gij {Pf ^, As) can be obtained upon taking the limit rrij —)■ 0 from the result of the previous 
section, 


f As 


mt 


( As 


G.,{Pt„As) = C{Pl„rr4) " ) iM 1 " 21 n( 7 ^ ) - ^ In | 


I 

\ in 


mj 


+ e{pf^-mi) 


2Li2 




(4.36) 


with the function C from Eq. (3.24). Similarly to the case of non-vanishing ruj, Gij{P'^, As) does not 
comprise an integral over the energy flow in the collinear regions, so that the quantity Qij (p 2 -,As,z) is 
not needed either. 

(Hi) Light emitter i and massive spectator j 

For a light emitter with mass we have to consider the hierarchy 


< As <C m|, P? . 


(4.37) 


Since the two limits m) <C As and As <C mf do not commute, the integral Gij{Pfj, As) does not result 
from the simple limit Wj —> 0 applied to the previous cases. An explicit integration delivers the result 


G.,iPlAs)=cl^,w^] - In 


+ 0 ( 7^2 



(4.38) 


It contains a term In^ As that results from a combination of the soft and the Tcollinear singularities, 
which are both isolated by the cut parameter As in the integration over the hard photon emission which 
is based on matrix elements with Wj = 0. 

In order to calculate Qij{Pfj, As, z) for non-collinear-safe observables, we first derive the functions 
2 / 1,2 (-z) for the integration boundary of the i/ij-integration in Eq. (4.28). The relevant range in yij is 
characterized by 0{m)/Pfj) < yij < 0{Ay), leading to the region defined by 


, , m) 1 — z 

= -d 


y 2 {z) = Ay. 


(4.39) 


The function Qij{P'lj, As, z) for the non-collinear-safe contribution is easily calculated to 


z) = ^ ^ In 

L — Z 


Asz 

m){l-z) 


2z 

1-z' 


(4.40) 
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Figure 4: Generic diagrams with a final-state emitter i and an initial-state spectator a (left), and with an 
initial-state emitter a and a final-state spectator i (right). 


fiv) Light emitter i and light speetator j 

In this case As) can be calculated upon taking the massless limit mj —)■ 0 of the previous 


ij 

case. The final result is 

GijiP'tp^s) = C{Pf^,mi) - In' 


■2 A i„2 ( _ 3. I 1 _ I ^ 

PfJ 2 U§] 2 


2 ’ 


(4.41) 


with Qij{P^j,As,z) as given in Eq. (4.40). 

4-2.2 Final-state emitter and initial-state spectator, and viee versa 

In the cases of a final-state emitter i and an initial-state spectator a, and of an initial-state emitter 
a and a final-state spectator i (depicted on the left- and on the right-hand sides of Fig. respectively), 
again the variables and abbreviations of Ref. 16 are used. The integrals over the full soft-collinear regions 
are defined as 


gupL^s) = 
G,,{Pl,As) = 
where the variables 


'Xq 

nX-i 


dxjo, 

p4 

la 

r-ziGic,) 

dz 

Uia 

^iaRiais^ ia) 

•t Zi{Xia) 


p4 

la 

r-z-iGic,) 

dz 

Uai 


J Zx(Xia,) 


PaPi+Pak-Pik 


d-Zja gfa (Pi^Pa, k) [OiOa + OtOa + OiOa), 
9a^'^\pa,Pi, k) [OiOa + O.Oa -f OiOa), 


PaPi +Pak 

are confined to the boundary determined by the variables 


PaPi 

PaPi +Pak' 


(4.42) 

(4.43) 

(4.44) 


■^ 1 , 2 ( 2 ;) = 


PfalPia - ^{Pia + 2^^)] =F ^P^i^ “ “ ^rn'^rrP^x'^ R,a{x) 


^PliPl - <Pfa - ^l)] 


2m jm-, 

a^i = 1 + p2 > 


(4.45) 


with xo being set by the least energy to trigger the process without photon emission]^ Here we use 

P^a=P^ + k-pa, Pf^ = Pf^-ml-mf-mi, X^a = (4-46) 


with A(a, 6 , c) defined in Eq. (3.10) and 




\/ {Pia + 2mlx)^ - 4Plmlx^ 


(4.47) 


^For > 0 and 0 < a /< rria — mi, the lower limit has to obey xq > x = 


'la. (^ rria - 


see Ref. 
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Note that we always have < 0, but Pf^ is not necessarily negative for massive particles i and a. The 
soft and collinear regions are characterized as follows, 


soft region: fc —>• 0 


i-collinear region: k Q, k \\ pi 


a-collinear region: k ^ 0, k \\ pa 


where 


Ax = 


As 



< As 


Xia > (1 + Ax) ^ 

< As 


^ia P 1 A 2 , 

< As 


Xia > (1 + Ax)“^ 

> As 


Zia <1- Az, 

> As 


Xia < (1 + Ax)“^ 

< As 


Zia>l- Az, 


(4.48a) 

(4.48b) 

(4.48c) 


Az = XiaAx. 


(4.49) 


By construction, the limits a;o,i respect the condition xq < 1/(1 + Ax) < Xi. Results for the integrals 
Gia{Pfa^ As) and Gai{Pfai ^s) for the various mass configurations of emitter-spectator pairs are listed in 
the following sections. 

Similar to the previous case of both emitter and spectator in the final state, the full integrals (4.421 
and (4.43) are not sufficient to deal with non-collinear-safe observables if the emitter particle f = i or 
/ = a is light, where <C As. To keep track of the momentum flow in the collinear region of the emitter, 
the Zia or the Xia integration has to be done numerically in the respective cases. To this end, we define 


Qia{Pfa^As,z) 

Qat{Pfa^As,x) 


fX2{z 


dx. 


Pt 


( 2 ) 2'\/ XiaRia (^ia) 


9tf°\pi^Pa, k) {6^9a + OiOa + OiOa), 


pi 


fZ2ix) 


k) {0i9a + 9i9a + 9,9a), 


2 \/^iaPia (^) Jz-i{x 


(4.50) 

(4.51) 


using X = Xia and ^ = Zia as shorthands. In the case of a final-state emitter, zpi is the fermion momentum 
after collinear photon emission from the final state; in the case of an initial-state emitter, xpa is the fermion 
momentum after collinear photon emission from the initial state, i.e. the momentum xpa enters the hard 
non-radiative process. The functions xi^ 2 (z) appearing in Eq. (4.50) result from interchanging the order 
in the integrations over Xia and Zia and are derived below for the relevant cases. The soft singularities 
in the integrals again are extracted from the numerical integration upon introducing plus distributions in 
the integrals. 


gUPL^s,z) = [gupL^s,z)]^ + 5 {i-z)G,a{PLAs), ( 4 . 52 ) 

ga,{Pl,As,x) = [ga,{Pl,As,x)]^+5{l-x)Ga,{Pl,As). (4.53) 

(i) Massive final-state emitter i and massive initial-state spectator a, and vice versa 

The case of massive initial-state particles is mostly interesting for decay processes, since the masses 
of colliding particles are typically much smaller than the scattering energies in high-energy collisions. If 
both final-state emitter i and initial-state spectator a (or vice versa) are massive, we have to respect the 
following hierarchy of parameters. 

As <^m‘f,ml,\P^J,\P^J. (4.54) 
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No collinear singularities are encountered, and the effective integration region in Gia,ai is a 

neighbourhood of the point {xia, Zia) = ( 1 , 1 ), so that the functions Gia{Piaj z) and GaiiPfa^ 
not required. The full integrals read 


GUPLM = ^ |ln(bi)ln t - 2 Li 2 ( 6 i) 

V L V / 


+ y-5ln>.) 




^ia 


("^) + 2 + ha API) T 0{PA ha API) for 

I in 1 nd n I + inicoOi j in | 
coj ^ 

A + A - An 


Li iif'ai 


,2 PL hJh\iJ”'y>' 




— In 


mi: 


ln( 6 i) — 2 Li 2 (co) + 2 Li 2 ( 5 i) + In(co) 


+ 2 In + farAPL) T 0{PA fa^APl) for ^ ma, 


with the auxiliary functions 


(4.55) 


(4.56) 
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2 mf 
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fai,±iPl) = ± 
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.Llnf In 


AA\ \2mA pItVaJ \-pItVA 
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p2 ' p2 


-2 Li, 


2 mi 


, Pi T 'J\a 


6 


Pi 


In 


2wf \ _ _ 2A , V^ l 

-pItVa) pi pi y 


(4.57) 


PI ± aa) 


(4.58) 


Here we partially used the results of Ref. 
/S.X —)• 0, and 


where we have identified (1 — xq) with Ax in the limit 


bi = 


2mf - Pi - 


P^ 


Co = 


AAa 


Pi - 


2mAPl + AAa’ 

(a) Light particle i and massive particle a 

If particle i is light and particle a is heavy, we have to consider the hierarchy 


(4.59) 


< As < TO„, |P,„|, |. 


|p 2 

I IC 


(4.60) 


so that we encounter a collinear singularity in the function GiaAl, but not in Gai{Pl, As). There¬ 
fore, we can obtain Gai{Pl, As) from Eq. (4.56) upon taking the limit rui —i- 0, but to obtain GiaAl, 
for small mi we have to perform a new integration. The explicit results are 


GAPlAs) = c 
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ml - Pf, 
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(4.61) 
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(4.62) 


20 












































In order to calculate the function QiaiPfa, z) for the non-collinear-safe contribution we have to perform 
the integration in Eq. (4.50) over the Xia range, which is bounded by 


1 — Aa; + < Xia < I -\- 




A simple integration yields 


Qia{Pfa, As, z) = In 

1 — z 


Pi. 


Asz 


+ a(mt/Pl). 


2z 


mf(l - z) 


1 - z 


(4.63) 


(4.64) 


finj Light particle a and massive particle i 

If particle a is light and particle i is heavy, we have to consider the hierarchy 

< As < ml, \Pfl\, \Pfl\, 

so that we encounter a collinear singularity in the function GaiiPl, As), but not in GiaiPfa, As). We can 
obtain Gia{Pl, As) from Eq. (4.55) upon taking the limit ma —t 0, but to obtain GaiiP^, As) for small 
TOa we have to perform another integration. The explicit results are 


(4.65) 
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(4.66) 


(4.67) 


In order to calculate the function Qai{Pl, As, x) for the non-collinear-safe contribution we have to perform 
the integration in Eq. (4.51) over the Zia range, which is bounded by 


I — xAx < Zia < 1 + 


A simply integration yields 


Gat{Pl, As, x) = LtH In 
1 — X 


As 


2x 


l{l-x)J 1-x' 


(4.68) 


(4.69) 


(iv) Light particles a and i 

Finally, we consider the case of both a and i being light, so that the mass hierarchy is 


< As < \P^J. 


(4.70) 


Here we always have Pi < 0. The full integrals Gia{Pl, As) and Gai{Pl, As) can be obtained from the 
results above where the respective emitter is light upon taking the limit of a massless spectator, i.e. upon 
taking —>• 0 in Eq. (4.61) and rui —>■ 0 in Eq. (4.67). The results are 

As 


G,a(Pl, As) = £(-Pl, ml) 
GalPl,As) = £(-Pl,ml)-hP 
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ta 
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-Pi 


-^In 

-5ln 
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-Pi 


TT^ 3 

~ Y ^ 2' 

7T^ 1 


(4.71) 

(4.72) 


The corresponding functions Gia{Pl, As, z) and Qai{Pl,As,x) for non-collinear-safe observables are the 
same as given in Eqs. (4.64) and (4.69), respectively. 
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Figure 5: Generic diagram with an initial-state emitter a and an initial-state spectator b. 

4-2.3 Initial-state emitter and spectator 

Although not needed for the treatment of decays, for completeness here we collect the formulas with 
both emitter a and spectator b in the initial state with masses rua and mb- The corresponding structural 
diagram is shown in Fig.[^ Using again the kinematical variables introduced in Ref. 16 , 


^ab — 


PaPb -Pak - Pbk 
PaPb 


Uab — 


Pak 

PaPb 


and 


s = s — 


^ab 


S = (Pa+Pbf, 

the integrals over the full soft-collinear region is defined as 

/■^l j. ,,2 pV2{Xab) _ _ _ 

Gab{s,As)= / dXab^r^J= / dyabgab \P‘>-^Pb,k){9aOb + 9aOb-\-0aOb) 

Jxq Aab Jyi{xab) 


(4.73) 


(4.74) 


(4.75) 


with the integration boundary defined by 
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_ 2m~^y/s 
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2 s 


{l-xf- 




(4.76) 


where again xq is the lowest kineamatically allowed value. The soft and collinear regions are characterized 
as follows, 


soft region: 

fc ^ 0 

1 

I Saj < As 
[Sh-y < As 


yab < Ay, 

^ab Uab 

Ay, 

(4.77a) 

a-collinear region: 

fc ^ 0, 

fc II Pa j 

I Saj < As 

1^ 5^7 As 


yab < Ay, 

^ab *^1 Uab 

Ay, 

(4.77b) 

6 -collinear region: 

fc ^ 0, 

k\\pb j 

I Sa7 > As 

1^ 5^7 ^ As 


yab > Ay, 

^ab Uab 

Ay, 

(4.77c) 


where 


Ay = —. 
s 


(4.78) 


Similar to the previous cases of emitter/spectator configurations, in addition to the full integral (4.75) we 
have to introduce the partial integral 


j.g2 |■V2(x) _ _ _ 

Qab(.S,As,x) = ^ ^ I dyabgabiPa,Pb,k){ 9 a 9 b + 9 a 9 b-\- 9 a 9 b) 


2-\/Aab . 


yiG) 


(4.79) 
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to keep track of the momentum flow in the collinear region of the emitter. Again we used x = Xab to 
shorten the notation. The momentum xpa is the fermion momentum after collinear photon emission from 
the initial state, i.e. the momentum xpa enters the hard non-radiative process. The soft singularities in 
the integrals again are extracted from the numerical integration upon introducing plus distributions in 
the integrals, 


Gabis, As, x) = [Gab{s, As, + 5(1 - x)Gab{s, As). (4.80) 

(i) Massive initial-state emitter a and spectator b 

If both initial-state emitter a and initial-state spectator b are massive, we have to respect the following 
hierarchy of parameters, 

< As < m^,m^,s, A (4-81) 

No collinear singularities are encountered, and the effective integration region in Gab{s, As) is a neigh¬ 
bourhood of the point (xab,yab) = (IjO), so that the function Gabis, As, x) is not required. The full 
integral reads 

" 7^ |ln(di) In j + 2 Li 2 (di) - y + ^ ln7rfi)| + 21n + 2 

fab,+ is) T 0i‘^rnl2ml - s) fab,±is) for rria ^ rrib, (4.82) 


with the auxiliary functions 
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s s 


(4.83) 


Here we partially used the results of Ref. 16 , where we have identified (1 — xq) with 2Ay in the limit 
Ay 0, and 

j, = »;+ypvy (4 84, 


s -f 2ml -t V '^ab 

There is no need to introduce the function Gabis, As, x) if the emitter particle is massive. 
(ii) Massive emitter a and light spectator b 
The hierarchy between the parameters is 


<C As <C m^, s, s, mb = 0. (4.85) 

Similar to the previous cases, the integral Gabis, As) can be obtained upon taking the limit mb —>■ 0 from 
the result of the previous section, 


A,) = £(., m^) - 21„ (^) In (7) - 2 In - 11„ (7 

+ <-(.-2™7i„77)+2l77)-^+21„(7 


2 - 


2 m? 


(4.86) 


Similarly to the case of non-vanishing m^ and mb, Gabis, As) does not comprise an integral over the 
energy flow in the collinear regions, so that the quantity Gabis, As, z) is not needed either. 
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(in) Light emitter a and massive speetator b 

If the emitter a is light and the spectator b is heavy, we have to consider the hierarchy 

ml As rril,s,s. (4.87) 

We encounter a collinear singularity in the function Gab{s, As), i.e. we have to perform another integration 
to obtain Gabis,As) for small ma- The explicit result is 


Gabis, As) = C - In^ 




+ e{s 




(4.88) 


In order to calculate the function Qab(s, As, x) we have to perform the integration 
Pab range, which is bounded by 


in Eq. (4.791 over the 


-^(1-x) < pab < Ay. 
s 

A simple integration yields 


Qab{s,As,x) 


1-x \ml{l-x)J 


2x 

I — X 


(4.89) 


(4.90) 


(iv) Light emitter a and spectator b 

Finally, we consider the case of both a and b being light, so that the mass hierarchy is 

ml «C ml, <C As s. 


(4.91) 


The full integral Gabis, As) can be obtained from the previous case upon taking the limit of a massless 
spectator, i.e. upon taking 0 in Eq. (4.88). The result is 


As\ 


As 


Gabis, As) =£(s, - In'" ( ^ j - ^ In ^ j - Y + ^. 


(4.92) 


The corresponding function Qabis, As,x) is the same as given in (4.90). 

5 Application to the top-quark decay 

As an application of the techniques described in Sects. and in the following we present the 
calculation of the QCD and EW radiative corrections to the top-quark decay at NLO. First, we summarize 
our calculational setup in Sect. 5.1 Then, we discuss the different contributions to the top-quark decay 
width in Sect. |5.2[ Finally, in Sect. |5.3| we present differential distributions in the top-quark rest frame 
for a set of kinematic observables. 


5.1 Calculational setup and input parameters 


Following Ref. 49 , the set of input parameters used in our numerical evaluation is given by 


mt = 173.34 GeV, 

= 80.385 GeV, 

= 91.1876 GeV, 
m^ = 105.6583715 MeV, 

G^ = 1.1663787 • lO"'^ GeV■^ 


mb = 4.78 GeV, 

= 2.085 GeV, 
rg® = 2.4952 GeV, 
Mh = 125.9 GeV, 
as(Mz) = 0.119, 


(5.1) 
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with the masses mt and Mh taken from Refs. 50 and 45 , respectively. Light fermions, i.e. all fermions 


but the top and bottom quarks, are considered massless. Further, we do not consider any mixing among the 
quark flavours in our calculation, i.e. we set the CKM matrix equal to the unit matrix. This approximation 
is justified because flavour transitions involving the third generation are practically forbidden (Vth ~ 1). 
Furthermore, any impact of a mixing among the first two generations cancels in the final result after 
taking the flavour sums owing to the unitarity of the CKM matrix and the approximation of vanishing 
quark masses for the first two generations. 

In order to describe the W-boson resonance, we employ the complex-mass scheme 


35 ■ 37 , a method 


that takes into account the effects of the instability and the off-shellness of the unstable particle in a gauge- 
invariant way, providing a consistent procedure at the one-loop level. The main idea of the complex-mass 
scheme is to consider the squared boson masses as complex quantities, defined through the gauge-invariant 
pole mass My and width Fy, 


iMvTv, V = W,Z, 


(5.2) 


and to introduce them directly at the level of the Lagrangian density. To apply the complex-mass scheme, 
the on-shell masses of the W and Z bosons given in Eq. (5.11 are converted to the corresponding pole 


masses according to 51-53 


My = 

Cv 


hy = 


IT 

Cv 


Cy = W 1 + 


/F^ 

[mOsJ 


V 


(5.3) 


To preserve gauge invariance, the complex masses have to be introduced in a consistent manner. In 
particular, this implies a complex definition of the weak mixing angle, 


= 1 - 


h-w 
Ai| ' 


and a generalization of the renormalization procedure for unstable particles, 
complex-mass scheme can be found in Refs. 


(5.4) 


Further details on the 


36 37 


The dependence on the bottom-quark mass is fully accounted for without approximation. The value 


of TOb given in Eq. (5.1) corresponds to the pole mass through the two-loop conversion formula given in 


Ref. 49 . Although all fermions with the exception of the top and bottom quarks are treated as massless in 


our calculation, a dependence on the mass of the light fermion arises when non-collinear-safe observables 


are considered (see Sect. 3.3). In the following, they are only considered for final-state muons where the 
experimental discrimination of a collinear muon-photon configuration is possible. Note that no event 
selection cuts are applied, so that the issue of collinear safety only arises in the event reconstruction. 
Thus, only differential distributions are concerned with this issue, while the integrated results, i.e. the 
(partial) decay widths, are collinear safe. 


The EW coupling constant a is derived from the muon decay constant via the relation 


aG„ = 


y/2M? 


w 


G,, 1 - 


Ml 


(5.5) 


This definition corresponds to the so-called G^-scheme, which avoids large logarithms induced by the 
running of a from the Thomson limit = 0) to the EW scale (Q^ ^ -^w) further absorbs 
dominant corrections to the p parameter into the coupling definition. Using instead of a in LO 
predictions changes the charge-renormalization constant by a term containing the quantity Ar 54 


which quantifies the radiative corrections to muon decay (see, e.g.. Ref. 55 for further details). Eor the 


strong coupling Og we employ the fixed value given in Eq. (5.1). 


5.1.1 Details on the calculations 

The implementation employing the slicing method of Section was performed using matrix elements 
that were calculated using the Weyl-van-der-Waerden spinor formalism (as formulated in Ref. 56 ) for the 


tree-level amplitudes and in-house Mathematica routines for the virtual one-loop corrections. Optimized 
phase-space parametrizations, based on Refs. 35 57 , were chosen for each of the pronunced structures 
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rLO[GeV] 

rNLOjQgY] 

(5“=[%] 

5“[%] 

semi-leptonic 

off-shell 

0.161065(1) 

0.148109(2) 

-9.38 

-Fl.34 


NWA 

0.163634 

0.150504 

-9.38 

-Fl.35 

hadronic 

off-shell 

0.483194(1) 

0.46242(3) 

-5.59 

-FI.29 


NWA 

0.490902 

0.46987 

-5.59 

-F1.31 

total 

off-shell 

1.449582(4) 

1.36918(6) 

-6.85 

-FI.30 


NWA 

1.472707 

1.39126 

-6.85 

-FI.32 


Table 1: Contributions to the decay width of the top-quark at NLO in as and a, divided into the semi- 
leptonic and hadronic decay channels. Results using the narrow-width approximation (NWA) are given 
for comparison. 


in the integrand induced by collinear poles or Breit-Wigner resonances. The numerical integration was 


performed using the Vegas algorithm 58 59 


A second independent calculation was performed using the subtraction method described in Section [3| 


Here, the matrix elements were obtained with the combined packages FeynArts 60 andFormCalc 6T|[^ 


An independent implementation of routines for the phase-space generation was employed. As in the 
calculation based on the slicing approach, the numerical integration was performed with Vegas. 

Both implementations resort to the Collier library 63 , which is mainly based on the results of 
Refs. 64,65 for the numerical evaluation of loop integrals and supports both dimensional and mass 


regularization for the treatment of IR singularities. 


5.2 NLO corrections to the top-quark decay width 

In Table [2 we collect the numerical results on the various contributions to the top-quark decay width, 
providing separate results on the semi-leptonic and hadronic decay channels. The relative correction 
factors S are defined with respect to the corresponding LO prediction T^®, 




pNLO QCD _pLO 


pNLO EW _ pLO 


(5.6) 


where pNLOQCD/EW ^.pg ]\jpQ prediction taking into account only QCD/EW corrections, while 

the full NLO QCD-FEW result is denoted T^*^. The QCD corrections to the semi-leptonic decay of the 
top-quark are known in semi-analytic form 
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^^ and are in full agreement with our numerical result]^ The 
full NLO predictior0 = 1.3693(2)GeV for the total width (see Table can be compared to the 

most recent CMS measurement [441 of 


Tt = 1.36 ± 0.02 (stat) (syst) GeV, 


(5.7) 


assuming m-t = 172.5 GeV, revealing good agreement between theory and experiment at the current 
accuracy level. 

A comparison between the results obtained with OGS and dipole subtraction is shown in Fig. As 
discussed in Sect. the slicing approach introduces a cut parameter As in order to isolate the singular 
region of the real-emission phase space. The dependence on this cut parameter cancels between the 
individual parts of the calculation in a non-trivial manner when combined to the final result, constituting a 
powerful check on the correctness of the calculation. Although this implies a certain degree of arbitrariness 


readjustment of the prefactor of the result given in Ref. 
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according to Gg —t G^ • Re(s^)/s^ is required in order 


to be consistent with the complex-mass scheme employed in our calculation. 


®The total width is calculated using lepton universality, i.e. the t mass m-r is neglected, since those effects are of 
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Figure 6: Dependence of the relative correction factor on the cut value As in the OCS method (red 
integration error bars) and comparison to the results obtained using the dipole subtraction method (band 
indicated by blue lines). The upper and lower borders correspond to one standard deviation of the 
numerical integration. The relative EW (left) and QCD (right) corrections are exemplarily shown for the 
case of a top-quark decaying semi-leptonically and hadronically, respectively. 


in the choice of the specific value for the cut parameter, in practice, sensible values for As are constrained 
to a certain range. This interval of validity is limited by the approximation used in the analytic integration 
on the upper, and by the stability of the numerical evaluation on the lower end. It is therefore necessary 
to perform a scan over the cut parameter to identify the region where both restrictions are simultaneously 
fulfilled and where the result can be trusted. Figure illustrates such a variation of As covering three 
orders of magnitude, clearly depicting the on-set of the aforementioned breakdown of the calculation. 
Furthermore, we observe good agreement between the results obtained by the two techniques which overlap 
within their Monte Carlo errors where the As variation develops a plateau. The OCS results presented 
in the following are obtained for a specific value of As where a point in the middle of the plateau was 
chosen after performing such a scan. This procedure is repeated for each decay mode (semi-leptonic and 
hadronic) and each type of correction (QCD and EW). 

Finally, we have compared the results of our calculation with off-shell W bosons to a narrow-width 
approximation (NWA) based on on-shell intermediate W bosons, which is defined by 


pNWA _ _ 
^ t^hff' - 


Tt—i-wb X 


BR 


w^//' 


Fw 


(5.8) 


at LO, where BR^_,. fp i s the branching ratio of the W boson. Note that we take the experimental result 
from Eqs. (5.1) and (5.3) as input for the total W width Fw, while the partial widths are calculated, in 
order to match the input procedure for the full top-quark decays with intermediate W bosons. At NLO, we 
linearize the corrections to in the corrections to the partial decay widths, in order not to spoil the 

proper comparison to the results for the full top-quark decays. Tableshows that the off-shell result and 
the NWA for the top decay widths differ by about 1.6%, in good agreement with the naive expectation 
of which quantifies typical off-shell effects of the W boson in integrated quantities. The 

relative corrections are influenced by the off-shell effects only at the sub per-mille level. Off-shell effects 
of the W boson will, however, play a more important role in differential distributions. 


5.3 Differential distributions in the top-quark rest frame 

In the following, we present differential distributions for the decay of the top-quark in its rest frame. 
After describing the event reconstruction in Sect. 5.3.H we present the numerical results for the semi- 
leptonic and hadronic decay channels in Sects. 5.3.2 and 5.3.3[ respectively. The differential distributions 
and their corrections are shown as plots that are subdivided into three frames with the following conven¬ 
tions: The upper panel depicts the absolute distributions, while the bottom two frames show the relative 
EW (middle panel) and QCD (bottom panel) corrections. 
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5.3.1 Event reconstruction 


It is not always possible to fully distinguish the particles of the process in the soft and collinear 
configurations. In our calculation, we try to be as close as possible to an experimental situation and 
distinguish the following scenarios: 


(i) Electroweak corrections to the semi-leptonic decay 

“Dresses leptons”: The EW corrections contain IR singularities originating from soft photons and 
configurations where the photon becomes collinear to the final-state charged lepton. For sufficiently 
inclusive observables these singularities cancel in the final result. To this end, we apply the photon- 
recombination procedure described in Ref. [ 66 | which attributes the momentum of a collinear photon 
to the charged lepton when close to it, i.e., when AR{£, 7 ) = Ay{l, 7 )^ -|- A(j){£, 7 )^ < 0.1, where 
Ay (A(/)) is the rapidity (azimuthal angle) difference between two objects. Such a recombination 
procedure is mandatory for a realistic experimental description in the case of electrons. 


“Bare muons”: In the case of muons, a collinear muon-photon pair can be experimentally resolved, so 
that no photon recombination is required. This results in non-collinear-safe observables, which are 


treated following the procedures described in Sections |3.3| and 4.2 


(ii) QCD corrections to the semi-leptonic and hadronic decays 


In contrast to the semi-leptonic case, where we assume that the bottom quark could be directly 
identified with a bottom jet, in the hadronic case we are confronted with the situation that experimentally 
indistinguishable partons are present in the final state. It is therefore necessary to define an IR-safe 
observable which groups soft and/or collinear partons into jets. In our calculation we adopt the inclusive 
generalized anti-fcx algorithm |67|g for two values of the clustering parameter, R = 0.6 and R = 1.0. The 


first value corresponds to the ATLAS default choice 69 , the second one will be motivated in Sect. 5.3.3 


We assume a bottom-tagging efficiency of 100% and identify the jet containing the b quark as the b-tagged 
jet. In the following, we refer to the jet candidate with the highest transverse momenturrj^as the leading 
jet. 


(Hi) Electroweak corrections to the hadronic decay 

When computing the EW corrections to the hadronic decay, we do not distinguish the photon from 
the other QCD partons and apply the same jet algorithm as described in the previous paragraph. 


5.3.2 Semi-leptonic decay of the top-quark 


The numerical results for the NLO corrections to the semi-leptonic decay of the top-quark, t —>■ b , 
are summarized in Fig.j^ They comprise differential distributions in the energy of the charged lepton (Eg) 
and of the bottom quark (Eb), the invariant mass of the bottom quark and the charged lepton (Mb^), 
and the cosine of the angle between the directions of the lepton and the bottom quark (cosOm). 

Following the reconstruction described in Sect. |5.3.1] we distinguish the two cases of “dressed leptons” 
and “bare muons” for the EW corrections. For observables that are not sensitive to the collinear £ ^ £"f 
splitting, i.e. which are collinear safe, both corrections must be nearly identical. This can be seen in the 
distributions in cos Oht and E-^ shown in the upper right and lower right plots in Fig. respectively. In 
the former case in particular, the photon is indeed recombined only if collinear to the positron, i.e. it does 
hardly change the direction of the charged lepton. The EW corrections have a minimum in the proximity 
of the peak (cos^bf ^ —0.8), are positive and smaller than 2%. For comparison, the QCD corrections 
reach down to —12%. The distribution in the b-jet energy is characterized by the fact that a 1 —)■ 2 
particle decay is effectively taking place when the W boson is on shell. In this case, the energy Ab of the 
b-jet tends to the fixed value 


e2^ = 


m( - -£- m 


2mt 


(5.9) 


®Note that no direction is distinguished in the top-quark rest frame and we have fixed an arbitrary reference axis to 
compute the transverse component. 
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Figure 7: EW and QCD corrections to the semi-leptonic top-quark decay as a function of Mb^ (top left), 
cos 0 b« (top right), Eg (bottom left), and (bottom right), for dressed leptons and bare muons. 


This explains the peak at Eb ~ 68 GeV. The asymmetric shape of the EW corrections is due to the 
typical radiative tail produced by photon radiation off a particle near a resonance energy. Most b quarks 
have an energy E^ « E^^ after the top-quark decay, so that photon radiation off b quarks, which is 
collinearly enhanced, reduces E^, leading to large positive corrections at energies below E^^. This effect 
is further enhanced in the relative corrections, because the LO distribution is strongly suppressed there. 
The same effect occurs in the QCD correction as well, leading to a similar shape, but even more than an 
order of magnitude larger. Note that in this case corrections more than 100% do not signal a breakdown of 
perturbation theory, they merely signal that the effective lowest order is provided by the radiative process 
with hard gluon (and partially photon) emission. A proper description of this tail, thus, would require 
the inclusion of multi-jet (and multi-photon) emission, as provided, e.g., by parton showers. 

For the remaining distributions, i.e. in Mbf and Eg, we observe EW corrections that are larger in 
magnitude for “bare muons”, stemming from the large mass-singular logarithm a In nif^ and reflecting the 
non-collinear-safe nature of the corresponding observable. The emitted photon in the £ —> 7 splitting 

removes some energy from the charged lepton. Contrary to the “bare muon” case, the “dressed lepton” 
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effectively can reabsorb some photons in its redefinition when the photon recombination procedure is 
applied. Hence, the EW correction to E^, always monotonically decreasing, are steeper for the “bare 
muon” than for the “dressed lepton”, as can be seen in the bottom left plot of Fig. For similar reasons, 
this effect is present also in the distribution in Mb^, as shown in the top left plot of Fig. Here, the 
EW corrections decrease from +6% to —4% for “dressed leptons”, while they drop from +8% to —10% 
when “bare muons” are considered in the top-quark decay. The QCD corrections to this distribution are 
decreasing down to —45% for large invariant masses. 

5.3.3 Hadronic decay of the top-quark 

Here, we present the results for the QCD and EW corrections to the top-quark hadronic decay width. 
In our setup, the corrections to the two processes t —)■ bud and t —>■ bcs coincide and are given in Table 
for the integrated decay width. 

In the top-left plot of Fig. |^we show the distribution in the invariant mass of the b-jet and the leading 
jet. In the upper frame we notice that the QCD corrections evaluated for R = 0.6 turn the distribution 
negative at Mbj ~ rrit. This is of course an unphysical effect: For a too low value of R, as in this case, the 
recombination between the real and the virtual parts is not sufficiently inclusive, so that the IR divergence 
leaves a trace as as InR corrections. This problem can be solved by enlarging the cone, e.g. to i? = 1, as 
shown in the same plot. With a wider cone, there is indeed a larger probability for the particles to be 
merged. In the case where all the final-state particles (excluding the b quark) are combined into a single 
jet, the invariant mass M^j reproduces the top-quark mass value. This explains the appearance of a peak 
at Mbj = rrif This obviously happens more frequently for a 1 —>■ 3 kinematics (as the LO) rather than 
at NLO. At NLO, the extra gluon can populate Mbj regions less probable at LO. This in turns explains 
the huge corrections in the lower frame of Fig. (top left): The magnitude of the NLO correction, shown 
in the middle frame, gets very large due to the strong suppression of the LO in the normalization. The 
EW correction is much smaller than the QCD correction, but is similar in shape. The curve for R = 0.6 
does not suffer from the same issue as the QCD counterpart. The R = 1 case is nonetheless shown for 
comparison. 

The distributions obtained choosing R = 0.6 or R = 1.0 are compared also for the other observables. 
The main difference is that the relative corrections induced by the choice of the larger radius are milder 
than those obtained when using R = 0.6, as expected from the inclusiveness argument. This effect can 
be clearly seen in the bottom-left and bottom-right plots in Fig. where the NLO corrections to the 
energy of the leading jet and of the b-jet are presented, respectively. The energy of the leading jet receives 
positive QCD corrections due to gluon radiation, which turn negative above the energy peak at ~ 65 GeV. 
As a consequence, the shape of the energy distribution at NLO is shifted to the left, more evidently for 
R = 0.6. In the low-energy tails the relative corrections are enhanced owing the normalization to the 
small LO prediction in this region. The QCD and EW corrections to the distribution in the energy of 
the b-jet, displayed in Fig. [^(bottom right), again show the distinctive radiative tails as already discussed 
for the semi-leptonic top decay in the previous section. 

The different choice of the jet-algorithm parameter affects also the the cosine of the angle between 
the b-jet and the leading jet, as displayed in the top-right plot in Fig. Once again, a larger jet radius 
implies a larger probability to merge all Hnal-state partons other than the b quark into a single jet. The 
latter will then be back-to-back to the b-jet since in this calculation the top-quark is considered at rest. 
This explains the peak appearing in the bin close to cos^bj = ~I when i? = 1 is employed. As in the 
previous case, the shape of the EW corrections match closely the QCD ones, however, are smaller in size 
by an order of magnitude. The strong rise of the QCD corrections toward larger values of cos 0bj again 
originates from the fact that gluon radiation populates the jet phase space more uniformly than at LO, 
where there is a strong tendency that the b quark and the hardest jet are back-to-back. 

6 Conclusions 

In this paper we have described extensions of two techniques for the treatment of infrared singularities 
in the computation of electroweak corrections at next-to-leading order: In Sect. we have supplemented 
the construction of universal subtraction terms within the dipole subtraction formalism to cover decay 
kinematics. The one-cutoff slicing method described in Sect. was derived on the basis of the dipole 
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Figure 8: EW and QCD corrections to the hadronic top-quark decay as a function of M\yj (top left), 
cos 9hj (top right), Ej (bottom left), and (bottom right), with j denoting the leading (non b-tagged) 
jet. The corrections are shown for the two choices R = 0.6 and R = 1 for the angular resolution parameter 
of the jet algorithm. 


subtraction formalism and organizes the singular contributions in terms of pairs of fermions. These 
results represent the first ones employing mass regularization and further allow for an arbitrary number 
of charged particles, a case that has not been discussed in detail in the literature so far. When dealing 
with more than three charged particles in the process, the appearance of new “surface” terms has been 
emphasized. For both methods we have covered all possible cases of massless and massive fermions in 
the final state and further worked out the respective extensions for the treatment of non-collinear-safe 
observables. The extension of the described techniques to the calculation of NLO QCD corrections is 
straightforward. 

Both methods have been applied to compute the NLO QCD and electroweak corrections to the top- 
quark decay width. In this computation, we have considered the b-quark as massive and evaluated the 
correction to the top-quark three-body decays fully accounting for the off-shellness of the intermediate 
W boson within the complex-mass scheme. Both the hadronic and semi-leptonic decay channels are 
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considered, further accounting for the case of non-collinear-safe observables in the latter. For the semi- 
leptonic decay, the NLO QCD corrections to the partial decay width amount to ^ —9.4% and ^ +1.3% for 
the electroweak corrections. In the hadronic decay channel, the NLO QCD correction are more moderate 
with ~ —5.6%, however, still much larger than the EW correction which amount to ~ +1.3%. While 
off-shell effects of the W boson modify the integrated decay widths by about 1.5%, the relative NLO QCD 
and EW corrections are modified only by few per mille by W off-shell effects. We have further presented 
differential distributions in the top-quark rest frame for different kinematic observables. 

The final result for the total decay width of the top-quark at NLO QCD+EW is given by = 

1.369 GeV, which corresponds to an overall correction with respect to the LO prediction of _ —5 5 % 
This breaks down to an overall NLO QCD correction of ~ — 6 . 8 %, and to a NLO EW correction of 
^ +1.3%, normalized to LO. As mentioned in the introduction, the total top decay width is known at 
NNLO QCD accuracy. The O^a^) correction, evaluated with a finite b-quark mass and a finite W-boson 
width in a soft-collinear effective theory and normalized to the LO, amounts to a further ^ —2.1% 27 29 . 
A comparison of our NLO QCD+EW prediction—with and without the additional NNLO correction of 
Refs. 


27 29 —reveals good agreement with the most recent measurement by the CMS experiment 44 


rCMS ^ ;l. 36 ± 0.02 (stat) (syst) GeV, 

FNLO ^ 3gg 


^NLO-l-NNLO QCD 


= 1.338 GeV. 


( 6 . 1 ) 


However, the quoted experimental error is still one order of magnitude larger than the typical EW correc¬ 
tions of ~ +1.3% addressed here. The experimental error is dominated by systematic uncertainties, such 
as the jet energy scale, where a significant improvement is difficult. The study of top-quark properties 
is an important part of the physics programme of future lepton colliders which aims for a precision at 
the order of 10—30 MeV in the experimental determination of the top width 70,71 and where the EW 
corrections are important. 

The codes developed in the context of this work can serve as a first step towards an unweighted Monte 
Garlo generator for the decay of the top-quark at NLO, which can be easily linked to any top-quark 
production process in order to simulate the decay at NLO accuracy. This is particularly important for 
processes with many particles in the final state when all intermediate particles are decayed. For instance, 
tt + V(V) {V = 7 , Z,W*), ttH, ttbb, and tttt, among others. Moreover, the top-quark decay width 
including both NLO QCD and EW corrections and W-boson off-shell effects comprises an ingredient in 
NLO predictions for the production of single top-quarks or tf pairs including top-quark decays and off- 
shell effects, such as the process pp —>■ bbW“''W“ —^ 6 /, for which NLO QCD corrections are known, but 
not yet the EW corrections. 

Appendix 

A Conversion from mass regularization to dimensional regularization 


In this appendix we briefly describe how to convert the results obtained using mass regularization 
(the default choice in the main part of this paper) to the corresponding expressions using dimensional 
regularization in D = 4 — 2e dimensions. To this end, we use the results of Refs. 18 47 for the general 


singular IR behaviour of one-loop amplitudes which also include the finite terms that arise from the non¬ 
commutativity of the limits e —>■ 0 and m —^ 0 and, thus, can be exploited for the conversion between the 
two regularization schemes. The following translation rules apply both to the functions G//'(P^, As) of 
the OCS method given in Section]^ and to the dipole subtraction functions and given in 

Section]^ and Refs. 16,19 , respectively. In the following we generically denote these functions Gffi. 

For an emitter particle with a non-vanishing mass m, i.e. m? ;+ As in the OCS method, only soft 
singularities arise which appear as lnm.y terms in mass regularization. The transition to dimensional 
regularization is simply accomplished by the substitution 


ln(' 




€ 




(A.l) 
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where we have defined the constant 


Q = (47r)T(l + e) = 


r(l - e) 


o 


(A.2) 


and fjL is the reference mass of dimensional regularization. 

If the emitter mass m is not light, but is used only as regulator in the limit m —> 0 {w? <C As in 
the OCS method), additional collinear singularities show up as Inm terms in mass regularization. The 
transition from mass to dimensional regularization is more conveniently described exploiting the fact that 
only the symmetrized sum G//' + Gf'f of the functions Gff appears in applications, since each charged 
particle can appear as emitter or spectator. In this symmetrized form the results of Refs. 18 47 can be 
exploited to derive the correspondence. However, we have checked that the following rule applies to each 
Gff separately as well. 

In mass regularization the IR divergences are always contained in the function C{P^,m?) defined in 


Eq. (3.24) with an appropriate kinematical invariant In dimensional regularization with a vanishing 


emitter mass m = 0 this singular function is universally given by 




1 3 , 


(A.3) 


Note that the difference between the results in the two regularization schemes, which is proportional to 
the combination £(P^, 0) — £(P^, m^), does not depend on the invariant P^, since the discontinuities arise 

from terms proportional to P’ which vanish in dimensional regulariation (to —> 0 

first), but not for mass regulzrization (e —>■ 0 first). 
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